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Abstract. Methods of calculation of loop diagrams in Heavy Quark Effective Theory 
(HQET) are presented. 



1 Hadrons with a Heavy Quark 

B meson is the hydrogen atom of Quantmii Chromodynamics (QCD), its sim- 
plest nontrivial object. In the leading approximation, the b quark in it just seats 
at rest at the origin and creates chromoelectric field. Light constituents (gluons, 
light quarks and antiquarks) move in this external field. Their motion is rela- 
tivistic; the number of gluons and light quark- antiquark pairs in this light cloud 
is undetermined and varying. Therefore, there are no reasons to expect that a 
nonrelativistic potential quark model describes B meson well. 

Similarly, Ab baryon can be called the helium atom of QCD. Unlike in atomic 
physics, where the hydrogen atom is much simpler than helium, B and Af, are 
equally difficult. Both have a light cloud with a variable number of relativistic 
particles. The size of this cloud is the confinement radius 1/v1qcd; its properties 
are determined by the large-distance nonperturbative QCD. 

The analogy with atomic physics can tell us a lot about hadrons with a 
heavy quark. The usual hydrogen and tritium have identical chemical properties, 
despite the fact that the tritium nucleus is 3 times heavier than the proton. Both 
nuclei create identical electric fields, and both stay at rest. Similarly, D and B 
mesons have identical "hadro-chemical" properties, despite the fact that b quark 
is 3 times heavier than c. 

The proton magnetic moment is of the order of the nuclear magneton e/(2mp), 
and is much smaller than the electron magnetic moment e/(2rne). Therefore, the 
energy difference between the states of the hydrogen atom with the total spins 
and 1 (hyperfine splitting) is small (of the order me /nip times the fine struc- 
ture). Similarly, the b quark chromomagnetic moment is proportional to 1/mb 
by dimensionality, and the hyperfine splitting between B and B* mesons is small 
(proportional to 1/mt,). 

In the limit m — > oo, the heavy quark spin does not interact with gluon field. 
Therefore, it may be rotated at will, without changing physics. Such rotations 
can transform B and B* into each other; they are degenerate and have identical 
properties in this limit. This heavy quark spin symmetry yields many useful 
relations among heavy-hadron form factors. Not only the orientation, but also 
the magnitude of the heavy quark spin is irrelevant in the infinite mass limit. 
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We can switch off the heavy quark spin, making it spinlcss, without affecting 
physics. 

The textbook [Tj contains a good introduction to QCD. Physics of hadrons 
with a heavy quark is discussed in the textbook . More examples of perturba- 
tive HQET calculations can be found in the lectures J3?4! ■ 

2 Heavy Quark Effective Theory 

Let's consider QCD with a heavy flavour Q with mass m and a number of light 
flavours. We shall be interested in problems with a single heavy quark staying 
approximately at rest. More exactly, let w <§; m be the characteristic momentum 
scale. We shall assume that the heavy quark has the momentum \p\ < lu and 
the energy \po — m\ ^ w; light quarks and gluons have momenta \ki\ < lo and 
energies |fcoi| < io. Heavy quark effective theory (HQET) is an effective field 
theory constructed to reproduce QCD results for such problems expanded up to 
some order in ui/m. 

In our class of problems, the lowest-energy state ("vacuum") consists of a 
single particle - the heavy quark at rest. Therefore, it is convenient to use the 
energy of this state m as a new zero level. This means that instead of the true 
energy po of the heavy quark (or any state containing this quark) we shall use 
the residual energy po = Pq — rn. The mass shell of the free heavy quark is 
Po = \J rn? + — m. At the leading order in 1/m, this means that po = does 
not depend on p. The free Lagrangian giving this dispersion law is Q+iZ?oQ, 
where Q is a 2-component spinor field describing the heavy quark at rest (we can 
also consider it as a 4-component spinor with the vanishing lower components: 
7oQ = Q)- Reintroducing the interaction with the gluon field by requirement of 
the gauge invariance, we arrive at the HQET Lagrangian 

L ^ Q+\DnQ + ■■■ (1) 

where all light-field parts (denoted by dots) are exactly the same as in QCD. 
The field theory Q is not Lorentz- invariant, because the heavy quark defines a 
selected frame - its rest frame. 

The Lagrangian Q gives the heavy quark propagator 

S{v) = ^:r^. S{x)^~m5{x). (2) 

Po + lO 

In the momentum space it depends only on po but not on p, because we have 
neglected the kinetic energy. Therefore, in the coordinate space the heavy quark 
does not move. The unit 2x2 matrix is assumed in the propagator It is 
often convenient to use it as a 4 x 4 matrix; in such a case, the projector 
onto the upper components is assumed. The static quark interacts only with Aq] 
the vertex is igv'^t'^, where — (1,0) is the heavy-quark 4-velocity. Loops of 
the heavy quark vanish, because it propagates only forward in time. 
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The Lagrangian |^ can be rewritten in covariant notations: 

Ly ^ Q^iv DQy H 

where the heavy quark field Qy is a 4-coniponent spinor obeying the relation 
"^Qv = Qy, and v'^ is the quark 4- velocity. The momentum p of the heavy quark 
(or any state containing it) is related to the residual momentum p by 

p — mv+p, Ip^l^m. 

The static quark propagator is 



S(P) 



1 + 



2 p-v + iO' 



and the vertex is igv^t"". 

One can watch how expressions for QCD diagrams tend to the corresponding 
HQET expressions in the limit m oo. The QCD heavy quark propagator is 



S(p) = — 2 ^r~^ — ^ -^^^ + o — 

p^ — m-^ Imp ■ V + p 2 p ■ V \m 



A vertex ig^^t"" sandwiched between two projectors may be replaced by 
\gvn''. 

Renormalization properties (anomalous dimensions etc.) of HQET differ from 
those of QCD. The ultraviolet behavior of an HQET diagram is determined by 
the region of loop momenta much larger than the characteristic momentum 
scale of the process w, but much less than the heavy quark mass m (which tends 
to infinity from the very beginning). It has nothing to do with the ultraviolet 
behavior of the corresponding QCD diagram with the heavy quark line, which 
is determined by the region of loop momenta much larger than to. 

The HQET Lagrangian (QJ possesses the SU (2) heavy quark spin symmetry. 
If there are Uh heavy-quark fields with the same velocity, it has the SU (2nh) 
spin-flavour symmetry. 

HQET has great advantages over QCD in lattice simulation of heavy quark 
problems. Indeed, the applicability conditions of the lattice approximation to 
problems with light hadrons require that the lattice spacing a is much less than 
the characteristic hadron size, and the total lattice length is much larger than 
this size. For simulation of QCD with a heavy quark, a must be much less 
than the heavy quark Compton wavelength 1/to. For b quark, this is technically 
impossible at present. The HQET Lagrangian does not involve the heavy quark 
mass TO, and the applicability conditions of the lattice approximation to HQET 
are the same in the case of light hadrons. 
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3 One-loop HQET Diagrams 

Now we shall calculate the one-loop HQET propagator diagram with arbitrary 
degrees of the denominators (Fig. 

r 

(-(fc+p)o -i0)"i(-fc2 -iO)"2 ■ 
It depends only on = Pa, and not on p. 

k 



k +p 

Fig. 1. One-loop propagator diagram in HQET 



The first step is to combine the denominators together. To this end, we write 
1/a" as 

1 1 /"^ 



a" rin) Jo 

(a- representation). Multiplying two such representations, we have 



r{ni)r{n2) 



-aiQi-a2Q2 „,"i-l„,"2-l 



an ^ a,^ dai da2 



We proceed to the new variables ai = ya, a2 = a, and obtain the HQET 
Feynman parametrization 



f r{ni + n2) y"'-^dy 



ai'aT r{ni)r{n2) Jo [aiy + 



^ ■ (4) 



■ a2\ 

If the denominator ai has dimensionality of energy, and 02 - of energy squared, 
then the Feynman parameter y has dimensionality of energy; it runs from to 
00. 

When combining the denominators in Q , it is more convenient to double the 
linear denominator; shifting the integration momentum fc — > fc — yv, we obtain 

^K + »2) f y"^-'dyd^k 

r{ni)r{n2) J [-fc2 + y(y _ 2uo) - 10]"^+"^ ' ^ ' 

The massive one-loop vacuum diagram (Fig. [21 is 



d-^k 



^d/2£M/2+^^^2)d/2- 



(-fc2 + - iO)" r{n) 
Substituting this integral with m? y{y — 2lu) into |(SJ), we obtain at < 

r{ni)r{n2) Jo ^ ^ ^ 
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We proceed to the dimensionless variable z = y/{—2u!). Then the substitution 
z + 1 ~ 1/x reduces this integral to the Euler S-function: 

. .d-n,-2n, rj-d + m + 2n2)r{d/2 - n^) 

^ r{-d/2 + m+n2) 

Our final result can be written as 
d-^fc 



Di^{k-v + uj)/uj, D2 = -k^, (6) 
r{-d + ni + 2n2)r{d/2 - n^) 



/(ni, 712) 



r{ni)r{n2) 



O 



Fig. 2. Massive one-loop vacuum diagram 

Ultraviolet divergences of this integral are the poles of r{—d/2 + ni + 2n2), 
and infrared ones - of r{d/2 — 712). The integral JHl has a cut at w > 0, where 
the real pair production is possible. If rii^2 are integer, /(ni,7i2) is proportional 
to /i, where 

/„ = ril + 2ne)r"(l - e) (7) 

is the combination of F functions which appears in the n-loop HQET sunset 
diagram. 

The one-loop propagator diagram of Fig. ^ in the coordinate space is just 
the product of two propagators. The HQET propagators in the p-space and the 
a:-space are related to each other by the Fourier transform: 

^" t-'eit), (8) 



- iO)" 2tt r{n) 

(-i)"+ir(n+l) 



e 







(-W - iO) 



n+l 



(9) 



For the massless ones, 



e-'^" d'p ^.^-2n-d/2 nd/2~n) 1 
(_p2 _ io)n {2^)^ r{n) (-a;2 + i0)<i/2-" ' ^ ' 



/2-n 

e'^-^ i^^- -..d-^n^dn ndl^-n) 1 



io)" r{n) (-p2 _ 



(11) 
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Multiplying the HQET propagator ||SJ| with the degree ni and the massless 
propagator ((Tn|l with the degree n2, we get 

r(ni)r(n2) 

(where —x'^ — (it)^)- Applying the inverse Fourier transform we reproduce 
the result ©■ 



4 Renormalization of HQET 

The HQET Lagrangian expressed via the bare (unrenormalized) quantities (de- 
noted by the subscript 0) is 

L = Qq'w ■ DqQo + ^ Qioi^lpo - ■mio)qiQ - ^Gq^^G'^'"' 

i 

where Do^Qq = (9^ — *5o^o^^")90: o-o is the gauge fixing parameter, and c is the 
ghost field. The rcnormalizcd quantities are related to the bare ones by 



1/2 

go = ZJ g , mio = Zmrrii , uq = Zau , 
where renormalization factors have the minimal structure 



(12) 



z-i + T^£ + (^ + T^jtS (13) 

All the renormalization constants in p2l) are the same as in QCD, where the 
heavy flavour Q is not counted in rif. In order to find the new constant Zq, we 
need to calculate the heavy quark propagator in HQET. 

The Lagrangian has dimensionality [L] = d, because the action S — J Ld'^x 
is exactly dimensionless in the space-time with any d. The gluon kinetic term has 
the structure (dAo)"^; hence, the dimensionality of the gluon field is [Aq] = 1 — e. 
Similarly, from the quark kinetic terms, the dimensionality of the quark fields is 
[<lio] = I ~ £• The covariant derivative Dq/x = dfj, — igaAg^t"' has dimensionality 
1, hence the dimensionality of the coupling constant is [go] — £• 

We define as to be exactly dimensionless: 



Here /i is the renormalization scale, and the factor exp(— 7^) is included for 
convenience. All the bare quantities, including go, are /i- independent. Differen- 
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tiating H14|) . we obtain the renormalization group equation for as{fJ-)'. 
d log as 



dlog/i 



-2e - 2P{as) , 

(15) 



where 



ldlogZ„ , ^ /as\ 

11 „ 4 



The bare (unrenormahzed) heavy quark propagator iS{uj) has the structure 
(Fig. 13 

iS{iLj) = iS'o(w) + iS'o(w)(-i)Z'(tj)i5o(ti') 

where Sq{lj) — l/w is the free HQET propagator, and the heavy-quark self- 
energy —1^7(0;) is the sum of one-particle-irreducible HQET self-energy diagrams 
(which cannot be separated into two disconnected parts by cutting a single heavy 
quark line). Summing this series, we obtain 

Siiv) = 1 . 



Fig. 3. Structure of diagrams for the heavy-quark propagator in HQET 

In the one- loop approximation (Fig. ^J, 

d'^k . ,, i . ,, — i / ^ k,,k 



After contraction over the indices, the second term in the brackets contains 
(fc • u)^ = {k-v + oj — Lu)^. This factor can be replaced by w^, because all integrals 
without k ■ V + Lu in the denominator are scale-free and hence vanish. Using the 
definition ©, we get 

^(^) = C^ ^%Iy/2'^ [2^(1' 1) + ^1 - 2)] , 



and, finally, 

y'{,.,\ - -n^^ , ,,,, _ , 

(4^)''/2 - 4 V d - 3 



^(-) - A ( «o - 1 - ) , (16) 
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Fig. 4. One-loop heavy-quark self-energy in HQET 

where Ii is defined in iQ. 

Therefore, with the one-loop accuracy, the heavy quark propagator in HQET 

is 



Now we re-express it via the renormalized quantities as(/i), a(/i)- Expanding 
the result in e up to the finite term, we get 

ujS[u) = 1 + Cf^^b-'^^' [i ^ a{^i) + Ae) where L = log ^ . 

47r£ /i 

This should be equal to ZQ{as{n), a{p))LijS,-{uj; fi), where the renormalization 
constant Zq^ag, a) has the minimal form (|13|l . and the renormalized propagator 
Sr{uj; n) is finite at e ^ 0. We obtain 

ZQK,a) - l + CF-^(3-a), (17) 

uS,{uj- = 1 -f Cf^^2 [{a{ii) - 3)L -f 2] . 
47r 

For any renormalization constant l|13|) , the corresponding anomalous dimen- 
sion is defined by 



dlogZ as , /"s\' 



It is more convenient to present results for anomalous dimensions instead of 
renormalization constants, because anomalous dimensions contain the same in- 
formation but are more compact. With the one-loop accuracy, it is sufficient to 
use the zeroth-order term —2e in 115|) when differentiating as(/i). We obtain 

7Q=2CF^(a-3) + --- 

Note this anomalous dimension vanishes in the Yennie gauge a = 3. 



5 Two-loop HQET Diagrams 

There are two generic topologies of two-loop HQET propagator diagrams, Fig.EjL, 
b. If one of the lines is shrunk into a point, the diagrams of Fig. [3;, d, e result. If 
any two adjacent lines are shrunk into a point, the diagram contains a no-scale 
vacuum tadpole, and hence vanishes. 
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1 2 
c 




1 2 
e 



Fig. 5. Two-loop HQET propagator diagram 



We write down the diagram of Fig. as 
d'^fci d-^fca 



n"i n"2 n"3 n"* n"^ 

1^2 ^3 ^4 ^5 



-Di = (fci +p) • w/w, i:>2 = (^2 +p) ■ w/w, 



(18) 



It is symraetric with respect to (1 <-> 3,2 <-> 4). If indices of any two adjacent 
lines are zero, the diagram contains a scale-free vacuum subdiagram, and hence 
vanishes. 

When one of the indices is zero, this diagram is trivial (Fig. If ~ 0, 
this is a product of two one- loop diagrams (Fig.jS]:): 



/(ni, n2, ris, n^, 0) I{ni,n-i)I{n2,ni) ■ 



(19) 



If m = (Fig. Eli), the (3,5) integral gives G(n3, n5)/(-fc|)"3+"=-''/2; this is 
combined with the denominator 4, and we obtain 



/(O, ^2,^3, "-4, "-5) = G(n3, n5)/(n2, r7,4 + ^3 + 715 - d/2) 



(20) 



(and similarly for 71,2 = 0). If 713 = (Fig. Et), the (1,5) integral © gives 
/(«!, ri5)/(— 2a;)"i+2"=^'*; this is combined with the denominator 2, and we ob- 
tain 

/(rii, 712, 0, 774, n-5) = I{ni,n<i)I{n2 -|- 77i -|- 2775 - d, 774) (21) 

(and similarly for 774 = 0). Here 0(771,772) is the massless one-loop propagator 
integral (Fig.EJ: 



d^fc 



i^''/2(V)''/'""'""'G(77i,772), 



G(7ll, 772) = 



{k+pf, D2 = -k\ 

r{-d/2 + ni+ ?72)r(rf/2 - ni)r{d/2 - 712) 
r{ni)r{n2)r{d - 77i - 772) 



10 



Andrey Grozin 




n2 

714 




14 + 713 + ns - d/2 



n2 

714 



X 



"1 "-2 ni n2+ni + 2n5 - d 

Fig. 6. Trivial two-loop diagrams 



71i 




712 



Fig. 7. Massless one-loop propagator diagram 



When all ti^ > 0, the problem does not immediately reduce to a repeated use 
of the one-loop formulae. We shall use a powerful method called integration by 
parts It is based on the simple observation that any integral of d/dki{- ■ ■ ) 
(or d/dk2{- ■ ■ )) vanishes. From this, we can obtain recurrence relations which 
involve /(ni, ^2, ns, 71.4, 715) with different sets of indices. Applying these relations 
in a carefully chosen order, we can reduce any /(tt-i, 71.2 , ns , 71.4 , 71.5 ) to trivial ones. 

The differential operator djdk^ applied to the integrand of H18|) acts as 



d 
9fe 



7l2 V 
D2 UJ 



"-5 , 



+ 1^2fc2 + ^2(fc2 - fci) 
D4 



(22) 



Applying {d/dk2) ■ (fc2 — fci) to the integrand of H18|) . we get a vanishing integral. 
On the other hand, from if^ . (^2 - fci) • v/uj — D2 ~ -Di, 2(fc2 — fci) • fc2 — 
-D3 — -D4 — -D5, we see that this differential operator is equivalent to inserting 



d — 7i2 — 'T.4 — 277,5 



under the integral sign [d comes from differentiating (^2 — fci)). Therefore, we 
obtain the recurrence relation |S] 



[d - 712 

where, for example. 



714 - 2715 + "22+1" + 7l44+(3" - 5")] / = , 



(23) 



1^/(711,712,713,714,715) = /(ill ± 1,112,113,714,715) . 
Expressing /(iii, 712, 713, 714, 715) from (|23|l : 

7i44+(5- - 3-) - 7122+1- 



/(7I1, 77,2,^3, 714,715) = 



d — 712 — 7lA — 27I5 



(24) 
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we see that the sum 711+713+715 reduces by 1. Therefore, applying sufficiently 
many times, we can reduce an arbitrary / integral with integer indices to a 
combination of integrals with 715 = (Fig. Et, CHl), iii = (Fig. Eli, (jSIB), 
713 = (Fig-Et, Thus, any integral /(tii, 712, 713, 714, 715) with integer rii can 

be expressed as a linear combination of and I2 , coefficients being rational 
functions of d. 

The second two-loop topology. Fig. Ebi 

"^'^^ "^'^^ -^^(-2^)2('^— J(71i, 712, n3, 714, 7.5) , 



Di ^ {ki +p) ■ v/uj , D2 = ik2 +p) ■ v/uj , D3 ^ {ki + k2 +p) ■ v/uj , 
I?4 - -kl , D5 - -fc| , (25) 

is much easier. If 714 = 0, or 715 = 0, or two adjacent heavy indices are zero, the 
diagram vanishes. If 713 = 0, this is a product of two one-loop diagrams (Fig.Et): 

J(7ii, 712,0,714,715) = /(ill, 714)7(712, 115) . 

If Til = this is a diagram of Fig. Et- 

j(o, 712,713,714,715) = /(713, 114)7(712 + 713 + 2714 - ^,"-5) ; 

the case 112 = is symmetric. The heavy denominators are linearly dependent: 

1 = D1+D2' D3 

When 711,2,3 are all positive, we can insert this into the integrand of H25f) . and 
obtain |H] 

J ^ {1- +2- -3-)J . (26) 

This reduces 711 +712 +713 by 1. Therefore, applying H26|l sufficiently many times, 
we can reduce an arbitrary J integral with integer indices to a combination of 
trivial integrals with 711 = 0, n2 = 0, 113 = 0. 

Using these methods, it is easy to calculate the two-loop heavy-quark self- 
energy (Fig. EJ). From this result, the two-loop rcnormalization of the heavy- 
quark field in HQET can be obtained 0: 

, , as /3a2 + 24a- 179 32 \/as\2 , , 

7Q = 2(a-3)Cp^ + Cp(^ CA + yrFnfj(^) +•.• (27) 

The method of calculation of three-loop propagator diagrams in HQET, 
based on integration by parts, has been constructed in [J] 



6 Heavy Electron Effective Theory 



Now we make a short digression into the abelian version of HQET - the Heavy 
Electron Effective Theory, an effective field theory of QED describing interaction 
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Fig. 8. Two-loop heavy-quark self-energy in HQET 



of a single electron with soft photons. It is obtained by setting Cf ^ 1, Ca 0, 
ffo — *■ Co, CKs ~* 0(. This theory was considered long ago, and is called the Bloch- 
Nordsieck model. 

Suppose we calculate the one-loop correction to the heavy electron propaga- 
tor in the coordinate space. Let's multiply this correction by itself (Fig. IHl). We 
get an integral in ti, t2, t'l, t'2 with Q < ti < t2 < t, Q < t'^ < t'2 < t. Ordering of 
primed and non-primed integration times can be arbitrary. The integration area 
is subdivided into 6 regions, corresponding to 6 diagrams in Fig.|51 This is twice 
the two-loop correction to the propagator. Continuing this drawing exercise, we 
see that the one-loop correction cubed is 3! times the three-loop correction, and 
so on. Therefore, the exact all-order propagator is the exponent of the one-loop 
term: 



S{t) = S'o(i)exp 



^(ii/2) 



2e 



(47r)'i/2 



-r{-e) ao - 1 



d-3 



(28) 



t'l 



ii 



t2 t 



t'. t 



■ + < 



■ + ■ 



• + 



Fig. 9. Exponentiation theorem 



In this theory, Za — 1, because there exist no loops which can be inserted 
into the photon propagator. Now we are going to show that = 1, too. To this 
end, let's consider the sum of all one-particle-irreducible vertex diagrams, not 
including the external leg propagators - the electron-photon proper vertex. It 
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has the same structure as the tree-level term: icqv^^F, F = 1 + yl, where A is the 
sum of all unrenormalized diagrams starting from one loop. Let's multiply the 
vertex by the incoming photon momentum g^. This product can be simplified 
by the Ward identity for the electron propagator fFig.lTTl: 

iSoip') ieow • q ^So{p) = icQ—^ [p ■ v - p ■ v] ^J— 

p' ■ V p ■ V 

= ieo So{p') - So{p) 

In the Figure, the gluon line with the black triangle at the end means that the 
vertex is contracted with the incoming gluon momentum q (it includes no gluon 
propagator!); a dot near an electron propagator means that its momentum is 
shifted by q. 



= eo 



Fig. 10. Ward identity for the free electron propagator 

Starting from each diagram for S, we can obtain a set of diagrams for A 
by inserting the external photon vertex into each electron propagator. After 
multiplying by q^, each diagram in this set becomes a difference. All terms 
cancel each other, except the extreme ones (Fig. Illll . and we obtain the Ward 
identity 

/l(cj,cj) = ; or 7 (CJ ,LJ) = ; . (29) 

The vertex function is thus also known to all orders. The charge renormaliza- 
tion constant Za is obtained from the requirement that the renormalized vertex 
function goFZ^^ Zq is finite. The factor Zq transforms in H29|) into 
and hence makes F finite. Therefore, the remaining factor (Z„Za)^/2 = 1 (this 
is also true in QED). In the Bloch-Nordsieck model, Za = 1 and hence Za = 1. 

Due to the absence of charge and photon field renormalization, we may re- 
place eo ^ e, flo — > a in the bare propagator H28|l . It is made finite by the minimal 
(in the sense of (|13|l 'l renormalization constant, which is just the exponent of the 
one-loop term 



Cq = exp 



(a -3) 



Aire 



and the anomalous dimension is exactly equal to the one-loop contribution 



7Q = 2(a - 3) 



47r 



Note that the electron propagator is finite to all orders in the Yennie gauge. 
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w///////m ^Mm^ 



+ 



eo 



• • • 



+ 



W////////M ^^^m 



W///Mm W/////////M 



eo 



J* 



W////////M W/////////M 



Fig. 11. Ward identity for the electron-photon vertex 



What information useful for the real HQET can be extracted from this simple 
abelian model? We can obtain the C| term in the two-loop propagator S{lS) 
without explicit calculation. There should be no C| term in the two-loop field 
anomalous dimension \H\ . 



Perturbative HQET 15 



7 Heavy-light Quark Currents 

Heavy-light quark currents are often useful (e.g., the weak 6 — > u current). QCD 
currents are expanded in 1/m; coefficients of these expansions contain HQET 
operators with the appropriate quantum numbers and dimensionalities. At the 
leading (zeroth) order in 1/m, we encounter the HQET heavy-light current: 

J(a*) = (asip)) JQ , Jo = qorQo , 

where -T is a Dirac matrix. 

Let the sum of one-particle-irreducible bare diagrams with the current ver- 
tex r, an incoming quark Q with residual energy w, and an outgoing quark q 
with momentum p, not including the external quark propagators, be the proper 
vertex r{uj,p) = F + A{iu,p) (here F in the right-hand side is the Dirac ma- 
trix). When it is expressed via the renormalized quantities as{fi), a{fi), it should 
become ZrF^{uj,p), where Zp is minumal Ijl^^f) . and the renormalized vertex 
Fr(u!,p) is finite in the limit s —> 0. When the proper vertex of the renormalized 
current r{ijj,p) is multiplied by the two external leg renormalization factors 

1/2 — 1/2 . ... / \ 1/2 ~ 

Zq Zq , it should give a finite matrix element. Therefore, Zj = ( ZgZq 1 Zp- 

Ultraviolet divergences of A{uj,p) don't depend on quark masses and the 
external momenta. Therefore, we may assume that all quarks are massless, and 
set Lo — 0, p = 0. An infrared cutoff is then necessary in order to avoid infrared 
1/e terms. In the one- loop approximation (Fis. I12|l . 



yl(0,0) = -iCF5o 



d'^fc j^jtv-" [gp. - (1 - ao)k^K/e 
d'^k 7o|^ - (1 - ao)fco 



{2ttY {k^fko 
Here ^ = kQjQ — fe • 7; the term with k yields after integration: 

1(0,0) ^ -iCpgia, ^ "^"^ ^ 



(2^Y (fc')' 



The ultraviolet 1/e pole of this integral is 
d'^fc 1 



{2ttY (fc2)2 



_ . 2^'^/2 
^^^\2^Ynd/2) 



00 

I,-l-2e 



dfc 



(47r)2e 1^ (47r) 



where i comes from the rotation to the Euclidean space, A is the infrared cutoff, 
and terms regular at e ^ are omitted. Therefore, 

Zr = l+aCF^. 

Aire 
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Taking into account Zq ()17|I and 



we see that the gauge dependence cancels: 

^ 2 47re 



Finally, 



1] 



47r 




Fig. 12. One-loop vertex 

We could equally easily obtain this result by calculating A[uj, 0) without in- 
frared cutoff. Keeping w 7^ is enough to ensure infrared convergence. At two 
loops this becomes mandatory. Calculating A{ijj, 0) in the two-loop approxima- 
tion (Fig. I13|l by the methods of Sect. |S1 we can obtain the anomalous dimen- 
sion [S] 



li = 



-— TT 

3 



(30) 



— TT 

3 
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Fig. 13. Two-loop vertex 



